Nucleation of breathers via stochastic resonance in nonhnear lattices 
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By applying a staggered driving force in a prototypical discrete model with a quartic nonlinearity, 
we demonstrate the spontaneous formation and destruction of discrete breathers with a selected 
frequency due to thermal fluctuations. The phenomenon exhibits the striking features of stochastic 
resonance (SR): a nonmonotonic behavior as noise is increased and breather generation under sub- 
threshold conditions. The corresponding peak is associated with a matching between the external 
driving frequency and the breather frequency at the average energy given by ambient temperature. 
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Intrinsic localized modes, often referred to as discrete 
breathers, have been the focus of a number of theoretical 
studies (for a recent review see e.g. [ij). One of the main 
reasons for the large interest in these modes is the fact 
that they provide a natural setup for energy localization, 
a paradigm of interest to many areas of physics as well as 
chemistry and biology. Such discrete breathers have been 
rigorously proven to exist as time-periodic localized exci- 
tations in nonlinear Hamiltonian lattices 0]. Ever since, 
they have been experimentally observed in a wide variety 
of different media, ranging from optical waveguides and 
photorefractive crystals to micromechanical cantilever ar- 
rays and Josephson junctions, as well as in Bose-Einstein 
condensates and layered antiferromagnets, among many 
others (EH. 

In the present work, having in mind realistic physical 
systems, we will focus on the nucleation of such breathers 
in a prototypical model system in which friction and am- 
bient noise are also present. This will lead us to the 
consideration of stochastic resonance (SR), the counter- 
intuitive phenomenon by which an appropriate dose of 
noise, instead of degradation, produces an enhancement 
of sensitivity of a nonlinear system to external forcing 
In spatially extended systems, SR has been shown 
for bistable Q or excitable media Q, usually associated 
with pattern formation 0], see e.g. the review of 

Our prototypical nonlinear model will consist of a 
quartic potential (the so-called hard lattice [§|) with 
the following equation of motion for each oscillator Xn 
{ji — 1, . . . , N) with mass m 

mXn = -U'{Xn) + k{Xn+l+Xn-l-2Xn)-aXn+S.n + Fnit), 

(1) 

where the overdot and the prime denote time and space 
derivative, respectively; k is the coupling parameter be- 
tween oscillators; a the damping constant; Fn(t) de- 
notes external driving; and ^„ is a Gaussian white noise 
with zero mean and autocorrelation {^nit)^m{t')) — 
2D5nm5{t — t'). The noise strength parameter D obeys 
the fluctuation-dissipation relation fc^T — D/a, where 
ks is the Boltzmann constant and T the ambient tem- 



perature. Thus, D is an indicator of the level of noise 
associated with thermal fluctuations. 

The on-site potential is given by U{x) — ax^ /2 + bx'^ / A 
with a,b > 0. This is a monostable nonlinear poten- 
tial which confers interesting properties to the lattice. 
The fact that is monostable implies that the noiseless 
lattice has a single stable steady state, allowing us to 
focus on the behavior of homoclinic orbits, rather than 
kink-like heretoclinic orbits, like in the bistable 4>'^ case 
(a < 0). This is one way in which our work differs 
from the standard bistable models of SR, or the inter- 
esting recent work of which considers noise-induced 
transitions between steady states in very short Fermi- 
Pasta-Ulam (FPU) type, boundary-driven chains. More- 
over, this model possesses a phonon band which can 
be easily determined upon linearization of ([1]) as G 
(v/f/"(0)/m, y/[U"{0) + 4fc]/m) ^. Thus, by exciting 
breathers with a frequency larger than the maximum fre- 
quency of the phonon spectrum, Wmax = 1/(0 -I- 4fc)/m, 
these localized modes are guaranteed not to interact 
with the phonons, avoiding effects like radiation reso- 
nances which slowly dissipate away the energy of discrete 
breathers [lol |. 

Indeed, the associated Hamiltonian system of this hard 
model (i.e., Eq. ([1]) without damping and noise) exhibits 
exponentially localized in space and periodic in time high 
(i.e., higher than Wmax) frequency excitations. An ex- 
ample is the dynamically stable Sievers-Takeno (ST) [ll| 
site-centered mode. As shown in Fig.[l] this mode is stag- 
gered, i.e., the neighboring sites oscillate out of phase. 
Hence, these breathers will necessitate that we impose 
staggered periodic driving in order to sustain them in the 
presence of dissipation. For this reason, we shall consider 
the following AC driving force 



F^it) - (-l)"^^osin(r!i). 



(2) 



Its effect on a deterministic FPU lattice has been studied 
in Refs. [13, 13 1. Furthermore, such staggered driving 



is straightforward to apply to driven nonlinear electrical 
lattices [2J|. Evidently, the driving frequency ft must be 



2 



1 


1 ' 


1 < u 


il , 1 



10 

n 



15 



20 



FIG. 1: Spatial profile of a staggered liigh-frequency (u = 3) 
Hamiltonian breatlier generated using metliods based on the 
anti-continuous limit 17]. The line is a guide to the eye. 
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chosen with the same value as the breather frequency we 
are interested to excite. 

The equations ([1]) are supplemented with (periodic) 
boundary conditions xq — and a;jv+i = 2:1- Gener- 
ally, <j)^ models has been commonly used to describe a 
wide range of phenomena from elementary particle colli- 
sions and cosmological domain walls IJ] to the behavior 
of an order parameter near the condensation point of a 
droplet in nucleation theory and phase transition phe- 
nomena [Tsj l. They also constitute a prototypical nonlin- 
ear model in statistical field theory [l6|. 

For the sake of a dimensionless description, we can 
scale all lengths and times by the characteristic length 
A — {a/by^^ and time r = {m/ay^^, i.e. i„ = a^n/A, 
i — t/r, which implies k — k/a, a = a/ {amy/"^ , 
D = bD/{a^my/^, and = Fn{b/a^y'^. In the fol- 
lowing we will omit the tilde symbols for simplicity. The 
dimensionless form of the equations of motion ([T]) then 
reads 

Xn — X^ X^ k{x^^\ -\- Xji—i 2x^) OiXyi 

Hn + {~irFosm{nt). (3) 

Let us discuss, for the moment, the deterministic be- 
havior of Q (i.e. without noise, D — 0). As mentioned 
above, the staggering factor (—1)" in ^ is important in 
order to sustain the breathers in the presence of dissipa- 
tion. When that factor is absent, the ST breather shown 
in Fig. [T] is always observed to die out, regardless of the 
driving strength. However, in addition to this staggering 
factor, the driving amplitude Fq must be large enough so 
that the energy supplied by the driving force is able to 
compensate the energy which is been dissipated by the 
friction force. Breathers in the dissipative deterministic 
lattice can be calculated using a method similar to the 
developed in • An analysis of those breathers shows 
that, for fc = 1 and a = 0.1 and a breather frequency 
of w = = 3, there is a threshold at Fth = 0.6092. 



FIG. 2: (Color online) Time evolution of the oscillator 
positions Un{t) ~ x„{t) -|- lOn for a system subject to a 
suprathreshold driving with Fo = 1 and f2 = 3, and noise 
strength D — 0.1 starting from a uniform initial condition 
Xn{0) = 0. Breathers are observed to form and vanish spon- 
taneously. The bottom panel shows a magnification of the 
time evolution of an oscillator which becomes the center of a 
high-frequency breather. 



When the driving amplitude is below this threshold (sub- 
threshold amplitude), the dissipation drains away grad- 
ually the system's energy and breather modes cannot be 
supported. In the following, unless explicitly stated, we 
will assume k — I, a = 0.1, and f2 = 3, which implies 

Wmax = \/5 < ^2. 

Having summarized the results of the deterministic 
case, let us now turn to the analysis of the effect of ther- 
mal fluctuations and how they modify this picture. Start- 
ing with the suprathreshold case {Fq > -Fth), we would 
expect that the addition of noise provokes the sponta- 
neous formation of breathers with the frequency dictated 
by the external driving, even if the system starts from an 
uniform configuration. This is indeed what is observed, 
as shown in Fig. [21 where a typical realization of the os- 
cillators' trajectories has been plotted. The frequencies 
and amplitudes of the observed noise-induced breathers 
are very similar to the deterministic ST breather depicted 
in Fig. [TJ The ST breather has the smallest size and 
the lowest energy of all possible collective modes at the 
frequency of the driver, thus being the most probable 
breather to be induced by noise. However, unlike the de- 
terministic case, the breathers are not permanent, and 
the same random force which has generated them is able 
to destroy them after a while. 

In the subthreshold case the driving force is not able 
by itself to sustain the breathers. Remarkably, however, 
the simulations show that these breathers are still pro- 
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FIG. 3: (Color online) Same as in Fig.[2]but for a subthresh- 
old amplitude Fo — 0.6 and noise strength D = 0.3. The 
bottom panel shows a magnification of part of the trajectory 
Xn{t) of one oscillator (boxed in the upper panel), as well as 
the time evolution of the center-oscillator of a deterministic 
{D = 0) and dissipative breather {a = 0.1) with suprathresh- 
old _Fo = 1 and same driving frequency (dashed-red line). 

duced in the system by the concerted action of noise 
and the driving force. This is shown in Fig. [3l The 
bottom panel shows a comparison of one of these noise- 
induced breathers with the deterministic stable breather 
obtained with a suprathreshold driving force. This ef- 
fect resembles one of the main features of the earliest 
manifestations of SR: the noise-induced hopping events 
of a Brownian particle in a bistable potential subject to 
a externally applied force which taken alone is not suf- 
ficient to produce the hopping events [4.]. Noise appears 
then as beneficial, enabling system transitions between 
the uniform background state and the coherent, local- 
ized breather patterns. On the other hand, we would ex- 
pect that if the thermal fluctuations are very large, they 
will destroy any degree of order in the system, just like 
in the mentioned Brownian particle example, producing 
the non-monotonic behavior characteristic of the SR phe- 
nomenon. This scenario should offer a generic mechanism 
for the formation of breather-like excitations in realistic 
nonlinear lattices where fluctuations are present. 

In order to quantify these effects, we have computed 
numerically the Fourier coefficient of an arbitrary oscil- 
lator at the frequency of the external driver, 

q^^ fj"" dte-^^'^x.it))^, (4) 

where the notation (. . .) indicates an average over noise 
realizations and the subindex oo indicates the long time 
limit of the noise average, i.e., its value after waiting for 
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FIG. 4: (Color online) The Fourier amplitude |q|^, de- 
fined in Eq. Q, as a function of the noise strength D for 
a suprathreshold driving Fq = 1.5 (solid line), a subthresh- 
old driving Fo = 0.6 (dashed line), and for a system with 
-Fo = 0.6 but without interaction between the oscillators, i.e. 
fc = (dotted line). 



time long enough for the transients to die out. Every 
time the oscillator xi is part of a breather, or a collective 
attempt to form a breather mode at the frequency of the 
driver, will result in a larger coefficient q. The choice of 
the oscillator in ^ is irrelevant because of the periodic 
boundary conditions. The simulations were performed 
with = 40 and = 60 oscillators, the larger value 
carried out in order to check that the results did not 
depend on the system size. 

We present in Fig. [4] the values of the amplitude \q\^ 
as a function of the noise strength D for a suprathresh- 
old and subthreshold driving. The non-monotonic be- 
havior with noise, a recognized signature of the SR phe- 
nomenon, is very clear in both the suprathreshold and 
subthreshold situations. Nevertheless, let us note that 
a non-monotonic behavior is also obtained when we in- 
hibit breather formation by canceling the interaction be- 
tween the oscillators. The dotted line in Fig. 2] shows this 
feature for a system with fc = 0. This is a well-known 
phenomenon for any single-oscillator system subject to 
a nonlinear monostable potential, both in the under- 
damped [13] or overdamped ji^l regime. In the former, 
the peak is produced by a matching between the driv- 
ing frequency and the intrinsic frequency associated with 
the most probable excitation ener gy, which is a mono- 
tonic function of the noise strength [l9|. In fact, a simple 
classical mechanics analysis of the undriven Hamiltonian 
system (a, D, Fq — 0) with k = shows that a periodic 
solution with frequency = 3 corresponds to an energy 
per particle value of E/N = 36. Under thermal equi- 
librium (i^o = 0); we need to raise the temperature to a 
value such that the noise strength is D = 1.67 in order to 
obtain an average energy per particle of the same value. 
It can be seen in Fig. 2] that this value of D is in good 
agreement with the position of the maximum. A similar 



4 



0.6 



0.4 



0.2 



' 1 ' 1 

of ^ 


1 ' 


•-•-^ X 

W X 


- 




X 




\ ^ 


f F =0 6^^- 

if 






^^--^^\x 




uniform driving 

1 



0.5 



D 



1.5 



V-Fo as 



FIG. 5: (Color online) Spectral amplification r\=\q 
a function of noise. Filled circles depict the linear response re- 
sults, computed numerically using Eq. (I?]), and the solid line 
the simulation results for Fo = 0.05. The dashed line cor- 
responds to = 0.6. Crosses and diamonds correspond to 
a linear response calculation neglecting correlations between 
oscillators separated by more than 2 and 3 sites, respectively. 
Filled triangles depict the linear response prediction for a sys- 
tem with a uniform external driving = Fq sin(nt) (for 
all n) instead of ([2]). The dotted line is a guide to the eye. 



argument also holds when we switch on the interaction 
between the oscillators, though now collective excitations 
associated with lower energies are allowed, resulting in 
larger amplitudes with a peak at a lower value of D. The 
total energy associated with the breather presented in 
Fig. [l]is i? = 30. If we consider that roughly only I — 5 
oscillators are needed to form the ST breather, an aver- 
age energy per particle of EjN — 6 would be necessary 
in order to optimize breather formation. This value is 
obtained at thermal equilibrium when D = 0.7, which is 
a reasonable estimation of the peak position observed in 
Fig. 01 A better estimation, namely a peak at D = 0.47, 
is obtained if we choose I = 7 oscillators. The appropri- 
ate value of I is determined by the correlation length, as 
discussed below. 

For an arbitrarily small driving amplitude Fq, linear 
response theory 2M, predicts 



{xi{t))oo= dsx{s)f{t- s), 
Jo 



(5) 



where f{t) = Fq sin(r2i) and xit) is the response function. 
This function can be expressed in terms of equilibrium 
time correlation functions. 



XW = -^|(^lWM(0))eq, 



(6) 



where (. . .)cq denotes equilibrium averages in the absence 
of driving and M{t) = J2n=ii~^)"'^ri{t) is part of the 
perturbing Hamiltonian —Mf{t) associated with the ex- 
ternal driving. The calculation of the Fourier coefficient 
of (ini), together with the definitions Xr = /g x(0 cos(51t) 



and Xi = /o°° x(^) sin(ri<), yields the following spectral 
amplification i] — i^i^ / p2 



(7) 



Fig. [5] shows the linear response results obtained from 
simulations at thermal equilibrium. It is seen that a driv- 
ing amplitude of Fq — 0.05 is small enough to be within 
linear response theory. Due to the presence of noise and 
dissipation, the correlation between two separated oscil- 
lators decreases rapidly as their separation is increased. 
In fact. Fig. [5] shows that a good approximation is ob- 
tained if the time correlation functions {xi{t)xn{0))cq 
corresponding to oscillators separated by more than 2 
sites (|n — 1| > 2) are neglected in Eq. This approxi- 
mation is consistent with the above energy matching ar- 
gument with 1 — 5, and the dominant dynamical role of 
ST breathers, leading to the peak estimation at D = 0.7. 
Larger driving amplitudes are associated with larger cor- 
relation lengths, producing a peak at slightly lower values 
of _D. A detailed quantification of the dynamical role of 
multibreather excitations [1, [13, [3] would be a natural 
topic for future studies. 

Finally, the calculated correlation functions allow us 
to readily compute the system response when the driv- 
ing force acting on each oscillator is not staggered but 
spatially uniform, i.e. M{t) — J2n=i ■'^nit)- The dot- 
ted line in Fig. [5] shows that the spectral amplification is 
very similar to that when breather formation was strictly 
forbidden (the case fc = shown in Fig. [4]), lacking the 
stronger peak at lower noise levels. As expected from the 
deterministic analysis, this uniform driving force is very 
ineffective in boosting breather formation. 

Summary. We have demonstrated the spontaneous 
formation of breathers in a hard c/)^ lattice subject to 
a staggered driving force. The SR mechanism is differ- 
ent from the one observed in the paradigmatic Brow- 
nian particle in a bistable potential. Here, instead of 
a stochastic synchronization between noise-induced hop- 
ping events and the frequency of the external driving, the 
formation of noise-induced breathers with the selected 
driving frequency is optimized for a given level of fluctu- 
ations such that the average energy per particle matches 
the intrinsic energy of the breather mode. Further, this 
effect is shown to be suppressed if a spatially-uniform 
driving is used instead of a staggered one. We expect 
that the results reported in this paper may be useful for 
the generation and manipulation of breathers in many 
practical situations, and to be independent of dimension 
and of the specific nature of the nonlinearity, provided 
that the model is monostable. 
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